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Abstract 

LHC searches for new physics focus on combinations of hard physics objects. In this work 

we propose a qualitatively different soft signal for new physics at the LliC the "anomalous 
underlying event" . Every hard LHC event will be accompanied by a soft underlying event due 
to QCD and pile-up effects. Though it is often used for QCD and monte carlo studies, here we 
propose the incorporation of an underlying event analysis in some searches for new physics. An 
excess of anomalous underlying events may be a smoking-gun signal for particular new physics 
scenarios such as "quirks" or "hidden valleys" in which large amounts of energy may be emitted 
by a large multiplicity of soft particles. We discuss possible search strategies for such soft 
diffuse signals in the tracking system and calorimetry of the LHC experiments. We present a 
detailed study of the calorimetric signal in a concrete example, a simple quirk model motivated 
by folded supersymmetry. In these models the production and radiative decay of highly excited 
quirk bound states leads to an "antenna pattern" of soft unclustered energy. Using a dedicated 
simulation of a toy detector and a "CMB-like" multipole analysis we compare the signal to the 
expected backgrounds. 



1 A New Signal: the Anomalous Underlying Event 



The LHC is about to probe the physics of the TeV scale. It holds great hope to discover new 
physics beyond the standard model and shed light on electroweak symmetry breaking and dark 

matter. The LHC detectors are designed to accurately detect several "physics objects" such as 
hard leptons, photons, jets, all of which are prompt, in addition to missing transverse energy and 
so on. Searching for new physics typically involves looking for an excess of events with a particular 
combination of such physics objects after an appropriate set of cuts is employed. Some of the 
conceptually simple searches include looking for an excess of events with leptons and/or jets and 
missing energy in models which include a WIMP dark matter candidate, or resonance peaks in 
the kinematic distributions of leptons in models with a new Z' gauge boson. The standard list of 
physics objects is also used for triggering. 

In addition to the standard set of physics objects, one can think of models which give rise 
to very non-standard objects. What is the value in these unconventional signals? Consider, for 
example, highly displaced vertices which were highlighted as potential signals for Hidden Valley 
models [1, 2, 3]. Discovery of highly displaced vertices can teach us that a new neighbor sector 
exists alongside the standard model. Studying the details of these events may teach us about the 
properties of the new sector. Another example is that of a stopped stable gluino that decays out of 
time in the hadronic calorimeter [4] . Such events can allow us to measure the lifetime of the gluino 
in split SUSY [5] and infer valuable knowledge about high energy parameters such a the SUSY 
breaking scale. Even though such signals are not what the detectors were designed to discover, 
non-standard analyses have shown that they may be discovered with a reasonable efficiency. If 
discovered, these signals may give us new and complementary information about the new physics 
at the TeV scale. 

In this paper we propose a new observable physics object at colliders - an anomalous underlying 
event. The underlying event (UE) is the often overlooked part of every LHC event. It is defined as 
everything in the event which is not the outgoing hard jets or leptons. Because of its omnipresence, 
the UE will be carefully studied and characterized for a variety of hard final states in terms of 
distributions of soft tracks and unclustered energy. Characterizing the UE will be a valuable tool 
for tuning Monte Carlo event generators. An anom,alous underlying event is one which contains 
a distinctively uncharacteristic distribution of soft tracks and diffuse energy when compared to a 
characteristic UE with a similar hard final state. Like any candidate signal, anomalous underlying 
events may and will be produced by SM backgrounds. Discovery of new physics with anomalous 
UE's requires a statistically significant number of anomalous UE's when compared to the SM 
expectation. Claiming such a discovery with confidence will require extracting the characteristics 
and of the typical UE and their variations from other LHC data. 

The signal and its search strategy are best explained in a concrete example. We will thus briefiy 
explain the model and signal which are the focus of this paper and sketch the search strategy for this 
example. We will defer both details and general lessons to subsequent sections. In this paper we will 
consider a new physics scenario in which some new particles with a mass of several hundred GeV 
are pair produced at the LHC. The new particles are charged under a new strong force and are thus 
confined to a single highly excited mesonic bound state. Due to the "quirky" [11] nature of the new 
dynamics the excitation energy may easily exceed several hundreds of GeV. After production, the 



2 




Figure 1: A schematic cartoon of the initial and final states of an LHC event with squirk production 
via an s-channel W^. The two protons are incoming along the horizontal axis. The squirks are 
produced and oscillate along the dashed axis. The final state includes an antenna pattern of 
soft photons (two cone like shapes aligned with the squirk production axis) and a pair of hard 
annihilation products, in this case. The search strategy will first involve discovering a resonance 
in Wj and then searching for signals of patterns of soft photons in the candidate signal events. 

excited bound state will emit soft radiation, and decay to the ground state, emitting many quanta. 
Some of these quanta will be soft photons which are emitted in a particular angular distribution. 
The ground state will then annihilate into a hard final state, for example a hard and a hard 
photon. The invariant mass of the VF+photon system reconstructs to the mass of the ground state 
meson (again, at several hundred GeV). All of the processes discussed above are prompt on collider 
time scales. A cartoon initial and final states of these events are depicted in Figure 1. 

The goal of the LHC search for this model would be to first establish that new physics is 
seen using standard hard physics objects emitted in the hard annihilation, and then to extract 
information about the nature of the new physics. In particular, detection of the unusual "antenna 
pattern" of soft photons in addition to the hard resonance will be a smoking gun signal of the 
strong dynamics and the presence of a bound state. What is a possible strategy to making these 
discoveries? In this case the existence of new physics may be demonstrated by a standard hard 
search. However, the correlation of new physics events with anomalous underlying event may 
teach us about the nature of the new physics, and perhaps enhance the confidence in the original 
discovery. A rough sketch of a search is shown in Figure 2 and described below. 

1. Establish the existence of an excess of events from new physics. 

Signal events will pass triggers with high efficiency due to the hard photon and lepton/jets 
from the annihilation. A promising search is to look for a peak in the W+photon invariant 
mass (or rather transverse mass) for leptonic W decays [13] . Due to the clean final state and 
the mass peak a signal-to-background ratio of order 1 may be achieved^ . 

^The signal to background ratio of order 1 may be achieved even when the transverse mass peak is smeared due 
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Figure 2: A schematic flow chart of a strategy for searching for anomalous underlying events in our 
example of squirk production. 



2. Identification of signal and control samples. 

We would like to determine if the new physics discovered in step 1 is associated with an 
anomalous excess in soft particles. For example, a resonance in VF+photon can be interpreted 
as a fundamental heavy W which contributes no new soft physics. If the origin of the BSM 
excess were a W' , the underlying component of the signal events is strictly a long distance 
phenomena and would thus be similar to the underlying component in SM background events. 
In order to test whether a new anomalous component in the underlying event appears in 
correlation with the hard signal, we can compare the underlying part of events in the signal 

to additional missing energy from hidden glueballs [13]. As will be described later, one may assume glueball emission 
is suppressed, in which case, the smearing of the peak is reduced and the signal-to-background will improve. In this 
work we will pick S/B~ 1 but the results can easily be rescaled to other signal-to-background values. 
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region to the underlying part of events that are known to be background dominated. 

One possible background dominated sample of events may be collected from neighboring bins 
in the kinematic distributions which were used to extract the signal, such as the regions 

immediately above or below the resonant peak (see Figure 2). Undclying event studies have 
shown that the characteristics of UE's are largely independent of the hard px in the event [6, 7] . 
The characteristics of the underlying event can thus be interpolated into the signal region 
reliably. 

In addition, since a large contribution to the underlying event at LHC is not directly related 

to the primary hard interaction, sample background underlying events may be collected from 
other channels with similar kinematics. For example, the UE in Drell-Yan + fi— with a \/l 
that is similar to the discovered resonance, may be studied relatively easily and compared 
with the signal. 

3. Compsirison of underlying events. 

Signals of anomalous underlying events may be discovered either in the tracking or calorimetry 
of the LHC experiments. UE studies often use the number density and pT distributions of soft 
tracks in the central region (or transverse regions for events with jets), e.g. [6, 7, 8]. In the 
case where the soft new physics is mostly in photons, the distribution of diffuse unclustered 
energy in the central EM calorimeter may be more promising. In particular, a multipole 
decomposition of unclustered energy may be used to recognize antenna-like angular patterns 
in the underlying event. As we will discuss, even in the case of soft photons a tracking signal 
may be seen from photons that convert in the tracking system. In this work we will consider 
both tracking and calorimetry signals, and present a more detailed study of the later. 

A similar strategy may be applied to other LHC searches, particularly searches for di-lcpton, 
di-jet, or di- photon resonances potentially teaching us that the new resonances are in fact de- 
excited bound states. An UE analysis of this type may also be applied to searches for large 
missing transverse energy. In this case an anomalous underlying component which is correlated 
with missing transverse energy may be a signal of light hidden valley particles that are decaying to 
soft SM particles and are "faking" the missing energy signal. 

In the example above the hard search alone gave a clear new physics signal and the UE study 
gave complementary information. One can also consider cases where a hard search is not likely 
to produce a significant excess above SM background (e.g. in di-jet searches). As we will discuss 
briefly, in that case the anomalous underlying event may be used as a tool to reduce backgrounds 
and improve the confidence in the initial discovery. 

The method of discovering soft new physics will depend on the channel in which it is produced. 
In the example discussed above, where a large number of soft photons are emitted, discovering the 
anomalous UE with tracking may be difficult (but perhaps possible) and calorimetry may be more 
promising. In this work we will focus mostly on this case. On the other hand, in some new physics 
scenarios the new soft physics is hadronic (e.g. in the case of the "hadronic fireballs" associated 
with colored quirks [11]). As we will discuss, in this case both tracking and calorimetry may be 
useful. 

Having sketched the goals of this work we will consider the example above in more detail. 
In section 2 we will present a quirk model which may lead to a large number of soft photons in 
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LHC events. We will briefly consider the dynamics and assumptions that lead to such events and 
parameterize some of the uncertainties. The model is directly motivated by folded supersymmetry, 
a recently proposed solution to the hierarchy problem using "superpartners" that are colored under 
a new gauge group instead of SM color. In section 3 we will describe how soft particle, particularly 
photons, interact in the LHC detectors and how they may be discovered either by tracking or 
calorimetry. In section 4 we focus on the calorimetric signal and present a simulation of a dedicated 
toy detector as well as simulations of the signal and background. We propose a possible observable 
based on a multipole decomposition of unclustered energy in the central calorimeter to discriminate 
signal and background events. We flnally give two examples of analyses that may be done in searches 
for anomalous underlying events. In section 5 we present a brief discussion and conclude. 



2 Quirks - New Physics with a Soft Diffuse Signal 

Before presenting a particular model, we will briefly discuss the following question - what type of 
soft signals can potentially lead to an observable anomaly in the underlying event? Obviously, an 
observable signal of this type will only exist if some new dynamics is generating a large multiplicity 
of soft particles which are unclustered. Additionally, an observable signal must dominate over the 
characteristic UE for the particular hard final state in question. Because most of the standard 
UE comes from QCD, most of the soft energy is emitted in the forward direction. In particular, 
distributions of energy and track multiplicity are expected to be flat in pseudorapidity (though this 
will ultimately be measured). Therefore new dynamics which emits soft particles into the central 
regions of the LHC detectors are more likely to be observable due to the lower background there. 
Finally, emission of soft particles in a particular angular pattern can help distinguishing signal from 
background. The new physics model we will now discuss will have all of these features. 

2.1 (S)quirks and their Dynamics 

We will now present a brief review of quirks, a new type of particle which exhibit new and interesting 
dynamics. These dynamics have been known for a while [10], and have received recent attention 
in the context of the LHC as well as their name in [11]. More recent discusions are in [?]. Similar 
dynamics may also be present in Hidden valley scenarios [1]. 

Though quirky dynamics may be exotic, it can arise from a seemingly mundane extension of the 
SM. Consider a model with a new strong force, QCD'. For concreteness we will take QCD' to be 
an asymptotically free SU(N) gauge force with a QCD' scale A of order a few GeV. We introduce 
some matter flelds, collectively called g', which are charged under the new force and also under 
the standard model interactions. As opposed to regular QCD in which there are quarks, u and d, 
whose masses are bellow the QCD scale, we will consider a case where all matter charged under 
QCD' is more massive than A 

m5/»A, (1) 

where even a little hierarchy between these two scales will be sufficient. When this condition is 
satisfied, the particles that are charged under the new force have been called "quirks" [11] (or 
squirks, depending on their spin) for reasons that will become clear presently. In particular we 
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will consider two types of quirks - uncolored quirks, which in addition to being charged under 
QCD' carry only electroweak quantum numbers, and colored quirks which also carry SM color. 
Beyond these two categories one can make various choices, such as the spin of the particles, that 
will matter only in the details. The examples wc will consider are motivated by a specific model, 
folded supersymmctry [9]. Starting with the uncolored case, consider the scalar squirk q' with the 
following quantum numbers 





SU(N)qcd' 


SU(N)qcd SU(2)i U(l)y 




N 


1 2 1/6 



This squirk doublet is quite similar to the doublet squarks in supersymmetric extensions of the 
standard model, with the exception that it is charged under QCD' instead of QCD. In fact the 

similarity to squarks may be exploited to solve the hierarchy problem in folded supersymmetric 
models in which the superpartners of SM quarks are charged under a new QCD' (with N=3) instead 
of QCD. 

Our choice of the QCD' strong scale is also motivated by folded supersymmetric models in 
which the two QCD scales are related by a Z2 exchange symmetry^. The lightest hadron in the 
QCD' sector is the glueball, most probably of the scalar type. This glueball is stable in a pure 
QCD-like theory, and may decay slowly to SM particles by higher dimensional operators. For our 
choice of strong scale, glueball decays occur far outside of any particle detector. 

A second particle motivated by folded supersymmctry that we may consider is a colored fermion 
quirk with quantum numbers 





SU(N)qcd' 


SU(N)qcd SU(2)i U(l)y 


~9' 


N 


N 10 



In folded SUSY such particles are the symmetry partners of the QCD gluon. The dynamics of the 
colored quirk share some common features with those of the uncolored squirk. In this work we will 
focus on the later and comment on the colored case occasionally. 

In addition to the modification of the hadronic spectrum, our QCD' sector will exhibit quali- 
tatively different dynamics at colliders once squirks are produced [10, 1, 11, 13]. Squirks will be 
produced at the LHC via weak interactions by either a Drell-Yan process or gauge boson fusion. 
They will typically be semi relativistic upon production, with a 7 factor of order two or so. The two 
squirks will fly off back-to-back but once their separation approaches they will start feeling the 
effects of QCD' confinement. More specifically the squirks will develop a gluon field configuration 
which, as the distance between them grows, will turn into a QCD' flux-tube or a QCD' string. 
In normal QCD with light matter the QCD string is torn promptly by producing a light quark- 
antiquark pair which can then separate the string to two (or many) pieces. Such a hadronization 
mechanism is possible because the energy density in the QCD string, is greater than the energy 
density needed to pair create, of order m^. In quirky QCD such a soft hadronization mechanism 
is obviously absent (i.e. exponentially suppressed [10]) due to the heavy quirk mass, A^ ^ m^,. 

^The different spectrum in the two QCD sectors will introduce a small logarithmic difference between the two 
strong scales through running. 



7 



Furthermore QCD' dynamics will not cause the heavy quirks to loose a significant amount of their 
kinetic energy, say by showering or fragmentation, over a distance of [21]. Instead the two 
heavy ends of the string will continue to move apart transferring kinetic energy into the QCD' 
string whose tension is a constant of order A^. The squirks will eventually stop when all of their 
kinetic energy was transfered to the string 

A^Lmax ~ Ek , (2) 

where Lmax is the maximum string length and Ej^ = s/i — 2mqr ~ m^/ is the initial kinetic energy 
in the squirk system upon production. The squirks will then be pulled together by the string, 
beginning oscillatory motion. We may thus conclude that squirk production at a collider is in fact 
a production process of a single highly excited mesonic bound state, "squirkonium" . As we will 
discuss in further detail in the next subsection, this excited meson will decay radiatively to the 
ground state and finally annihilate to hard decay products. For our choice of A all of the processes 
above will be prompt on collider time scales. 

What is our signal? The soft radiation emitted during the decay from the highly excited to the 
ground state is the soft signal which we will ultimately detect. In the case of uncolored squirks the 
signal will consist of many unclustered soft photons^. In the colored case the soft signal (which 
was dubbed a "hadronic fireball" [11]) will mostly consist of a large multiplicity of pions. The hard 
annihilation products will provide the hard primary signal. 

What are the hard primary signals to search for? In the colored case: the hard signal was 
discussed in [11]. The most common visible annihilation is to jets, but given the higher production 
rate for colored quirks, leptons or photons may be a feasible final state. 

For uncolored squirks: a neutral squirkonium state will typically decay to two hard (but 
invisible) QCD' glueballs which will not give a triggerable signal. However, the dominant production 
of squirks will be via an s-channel and thus the produced squirkonium is charged under standard 
model QED. Its decay products will then always contain a charged particle and leave a visible 
signal in detectors. An interesting decay channel which may dominate is W'^'y [13]. Because the 
annihilation will typically occur at or near the ground state, a resonance peak is expected in the 
invariant mass of Wj. 

The details of the search for the hard annihilation of uncolored squirks is discussed in [13]. 
Here wc will ask the following question: Civen an excess of events with a signal-to-background 
ratio of order one, can we learn more about the nature of new physics by carefully analyzing the 
underlying event? To answer this we will discuss the energy and angular distributions of the soft 
photon radiation in further detail. 

2.2 Soft Radiation 

We will now consider the dynamics of a highly excited squirkonium bound state and its decay. 

Because the bound state is produced with a very high principal quantum number, the system 
may be treated semi-classically [11]. The system can then be modeled by two heavy particles of 

^We will discuss soft hidden glueball radiation in the next subsection. 
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masses mi^2 = with an attractive linear potential V = A^jxi — X2\- In the center of mass 
frame the motion is parameterized by the total energy E and by the angular momentum I. The 
energy is initially set by the energy of the collision and will decrease as the oscillating system 
radiates. The angular momentum number / is initially of order one, and will change as the system 
radiates by one unit per radiated quantum on average. The angular momentum will therefore 
increase on average by a random walk. Classically, the angular momentum is related to the impact 
parameter h ^ 1/ ruq' , and is important to determine the likelihood of short distance effects such as 
re-annihilation as will be explained below. 

The classical trajectory of the two particles is easily computed and is to a good approximation 
linear oscillatory motion. The period of the motion is 



Pnear n^-q' 



P 



A2 



where Pncar is the the momentum of the quirk at the point of nearest approach (or at production 
for the first periods) and p ~ is the force acting on the quirks. 

The two particles arc charged under standard model QED (with charges ±2/3 and ±1/3) and 
are also triplets under QCD'. The accelerating charged particles will radiate to both of these sectors. 
As a simplifying first step we will completely ignore radiation of QCD' glueballs, returning to this 
assumption later. In this simple approximation the excited system will slowly decay by classically 
emitting electromagnetic radiation in accordance with Larmor's law 

E= :^p^. (4) 

The spectrum and angular distribution of photons may also be readily calculated by Fourier 
decomposing the retarded potential far away from the source [22]. Here we will discuss the 
qualitative features of the distribution, leaving details for Appendix A. The Fourier series will 
be dominated by the frequency of oscillation of the quirk system which for most of the motion^ is 

CO ~ (5) 

For a squirk mass of 500 GeV and a QCD' string tension of (5 GeV)^ the photon spectrum is 
dominated by frequencies of order a few hundred MeV. Such a soft photon by itself is not observable 
at the LHC. however as we will show, if a significant amount of the squirkonium's energy is lost 
though such photons, they may give an observable modification to the underlying event. 

The angular distribution of the emitted radiation is particularly important. A relativistic 
charged particle which decelerates and is brought to rest is known to radiate in a cone in the 
forward direction. A particle accelerating from rest radiates in a similar pattern. The oscillatory 
motion may be thought of as a succession of such accelerations an decelerations. The "antenna 
pattern of radiation" (see figure 1) may be understood as a combination of two such cones back to 
back. The details of computing this distribution are given in Appendix A. As we shall see later, a 
calorimeter surrounding the decaying system will thus see two circles of radiation in the r]-(p plane. 
This peculiar pattern will be the smoking gun signal of such radiative decays. 



''As the system looses energy, pnear will decrease as sJE^ — vn? and the frequency will gradually increase, however 
most of the energy will be radiated at frequencies near the initial one, of order IS? /niqi. 
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Before demonstrating how to search for these patterns we should consider some concerns and 
caveats. For example in order for a large amount of energy to be emitted in radiation we must 
ensure that the squirks do not annihilate in a highly excited state. This, however, is unlikely. To 
argue this we will adopt the semi-classical formalism of [11]. 

The probability for annihilation in a single crossing depends strongly on the angular momentum 
of the bound system, peaking at low Vs. Classically it is clear that annihilation will only be likely 
at low impact parameter, where a naive expectation would be that the annihilation probability will 
scale like aann/b'^ oc A more rigorous partial wave analysis gives a much stronger dependence, 
scaling like /^^~'. Given that the angular momentum will grow on average as radiation is emitted, 
the likelihood of early annihilation is determined by a "horse-race" between the annihilation cross 
section and the radiation rate. Very naively one can expect that radiation would win since the 
radiation rate is proportional to a and the annihilation cross section scales like a^. However to 
make a clear determination a more careful estimate is required. 

The probability for emitting a photon per period is given by 

Pr^^T-r^a. (6) 

and thus a photon will be emitted on average once per crossings. This can be compared with 
the probability of annihilating at low Vs. In our case of charged squirkonium, it was shown [13] 
that at high relative velocities, right after production, annihilation goes dominantly to SM fermions 
with a cross section 
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a = NQCB'NfV.ei^ (7) 

where Nf = 12 is the number of SM final states, aw is the weak coupling constant and Vrei ~ 1 is 
the relative velocity. Following [11], the probability for annihilation per crossing at high velocities 



is5 

,2 



P 



at;.ei~ArQCD'iV^^(^)'~10-^ (8) 

where 2mq/E is typically of order 1/2. Therefore, an order of 10^ photons will typically be emitted 

before the initial probability to annihilate becomes appreciable. However, the change in angular 
momentum due to the emission of these photons (by Al ^ 10 on average) will significantly reduce 
the probability to annihilate. 

This argument implies the bound squirkonium system will likely decay radiatively to a low 
lying state before annihilating. This is fortunate not only because the soft radiation provides an 
interesting signal, but also because the hard annihilation will occur at a fixed invariant mass, easing 
its identification independently [13]. 

Another caveat we should discuss is the possibility of loosing energy by the emission of QCD' 
glueballs. In the case we are considering hidden glueballs decay outside the detector and do not 

lead to an ol)scrA'al)lc signal. In fact, one may argue that glueball emission is likely to dominate 



^ In [11] it was argued that EM radiation is not likely to prevent a fermionic uncolorcd squirkonium system from 
annihilating early. This is because in the case studied there Pann was only a factor of a few smaller than a. This, 
however, highly depends on the details of system at hand. In particular the annihilation cross section for our scalar 
system is significantly smaller. 
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over radiation of photons because of the strong couphng in the QCD' sector. However, because of 
the gap in the QCD' spectrum, glueball emission may be suppressed. In the semi-classical limit 
the oscilating suqirks would prefer to radiate at frequencies that are in tune with the oscillation 
frequency uj ^ h? /mqi. No hadron in the QCD spectrum is light enough to be radiated in this 
frequency and we thus expect photons to dominate over the majority of the oscillation time. When 
the quirks are close to one another the semi classical limit breaks down and quantum emission 
of a glueball is possible. However, the probability for such an emission may be kinematically 
suppressed as well. According to lattice studies the glueball mass is expected to be of order 3-4 
times A [20]. Given the small hierarchy between the glueball mass and the QCD scale one may 
estimate this emission perturbatively in which case the probability for emitting glue is down by 
(rrtgiue/A)^ ~ 10~^ [11]. We can thus naively estimate that the rate of energy loss due to glueball 
emission is of order a few times 10~^A per period. 

How does this compare with photons? Recall that the probability for emitting a photon of 
energy h? /mqi per period is of order a. Very roughly the competition between photons and glueballs 
reduces to comparing the kinematic suppression of 10~^ to aA/mg/. Using A = 5 GeV and mqi = 
500 GeV the later is of order 10~^. We thus estimate that losing 10% of the energy into soft 
photons is quite reasonable. This however is a very rough expectation, and various factors may 
affect the answer in cither way. For example, in the large limit glueball emission is known 
to be suppressed as A^^^, giving an order of magnitude suppression even for Nc = 3. Given the 
high level of uncertainty in the fraction of excitation energy lost to photons we will take a more 
phenomenological approach. We will consider two cases: case (1) in which glueball emission is 
highly suppressed and all of the energy is emitted in photons, and case (2) in which we will only 
allow 10% of the energy to be emitted in visible photons and the rest is lost to glueballs. We will 
asses the prospects of observing a deviation from the standard underlying event in both of these 
cases, finding interesting results in both. 

Finally, we will consider the case of colored quirkonium. The story here is qualitatively similar: 
quirks are pair produced (with a larger cross section) and form a highly excited bound state. When 
the ends of the string are colored energy may be lost by radiating light QCD pions. Any kinematic 
suppression which may be present in the QCD' case will now be absent and we expect pions to 
dominate the energy and angular momentum loss. In [11] it was argued the produced bound state 
will ring down to a low lying state before annihilating in a significant fraction of production events. 
It may thus be feasible to search for the hard annihilation products of colored quirkonium states 
in association with a large multiplicity of soft pions. These new pions have been dubbed "hadronic 
fireballs" [11], but they may be thought of as a new additional component to the underlying event. 
As we will see, the techniques for searching for soft photons will also be useful for the hadronic 
case. With this in mind we are now ready to ask how a large number of soft particles, particularly 
photons, can be detected in an LHC detector. 



3 Soft Photons in an LHC Detector 

In this section we will explain how LHC experiments can detect and measure diffuse soft photons 
and which parts of the detector arc relevant for these measurements. Wc will find that a large 
number of soft photons can potentially leave observable signals in both the tracking system and 
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calorimeters of the LHC experiments. We will then consider the prospects and techniques for 
observing both tracking and calorimetric signals and distinguishing them from backgrounds. 



3.1 EM Showers and Soft Signals 

We will now briefly review how soft photons interact with matter triggering EM showers. A 

useful review of these subjects is in [14]. High energy photons (>50 MeV) interact with matter 
dominantly by conversion to an e^e~ pair. High energy electrons loose energy mostly by emitting 
bremsstrahlung photons. The successive repetition of these processes is an electromagnetic shower. 
Because photon conversion and bremsstrahlung emission are related by a crossing symmetry the 
shower may be characterized by one length scale, a radiation length or Xq. The radiation length 
Xq, is the typical distance over which an electron looses an order one fraction of its energy by 
bremsstrahlung. It is also 7/9 times the mean distance a photon travels before converting. Xq 
depends on the medium which the particles traverse, and is roughly 9.3 cm in solid silicon. 

As the EM shower progresses, the individual particles in the shower loose energy until a critical 
energy, Ec, is reached. Below E^. other energy loss mechanisms such as ionization dominate, quickly 
ending the shower. The critical energy is roughly Ec ~ (800 McV)/Z in material with an atomic 
number Z. Once a photon or electron have been produced below the critical energy they are 
promptly absorbed and that fraction of energy is considered "lost" or deposited. 

The soft photons we are considering will initiate EM showers in the LHC detectors. For our 
study we will need to know where inside the detector most of the energy is deposited and what 
fraction of the energy reaches the calorimeter. Given a shower that was initiated by a photon of 
energy Eq the fractional energy deposition is maximized after a distance 

Xmax = Xo (^In (^^^ + 0.5^ (9) 

The average fractional energy loss of an electromagnetic shower may be expressed as [14] 

liE^iHT^ (10) 

Eo dt T{a) ^ ' 

where t = xjX^ is the distance in units of radiation lengths. The parameters a and 6 are defined 
as (a — l)/6 = XmaxlX.Q and h ~ 0.5 for our purposes. This fractional energy deposition is plotted 
in Figure 3 for a shower initiated by a 0.1, 0.3 and 1 GeV photon in silicon. A rough material 
budget of the various components of the LHC experiments is also shown. It is interesting to note 
a slight difference between ATLAS and CMS. At ALTLAS the calorimeter is behind the coil and 
the cooling system, which dominate the pre-calorimeter material budget. At CMS the calorimeter 
is inside the coil. This difference is not crucial for photons in the energy range of interest. We can 
summarize the matter effects by the following rough estimates 

• Roughly 30% of the emitted photons undergo their first conversion to an e^e' pair inside the 

tracking system leaving to soft tracks. 
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CMS □ 




t = x/Xq 

Figure 3: The fractional amount of energy deposited by an electromagnetic shower as a function 
of distance traveled in radiation lengths for photon initiated shower with energies of 0.1, 0.3 and 1 
GeV. A rough estimate of the material budget of the two experiments in the central region (77 < 0.6) 
is shown. 

• An order one fraction (30-90%) of the energy reaches and will be deposited in the electro- 
magnetic calorimeter. 

This suggests that anomalous underlying events which are dominated by soft photons may be 
searched for, both in the tracking system and in the calorimeters. In particular, one can search for 
either or both of the following 

1. A high multiplicity of charge tracks in the central region. 

2. Soft energy deposition in the EM calorimeter with distinct geometrical distributions. 

Both of these searches should focus on the central regions where backgrounds are low and (for the 
calorimetric signal) the material budget before the calorimeter is minimal. We will briefly discuss 
the tracking signal in the next subsection before going into the calorimetric signal in more detail. 

Before we do so, we will consider anomalous underlying events which are dominated by SM 
hadrons, as we expect would be generated by colored quirks. In this case one would mostly expect 
charged and neutral soft pions. The neutral pions, which are roughly one third of the total, will 
promptly decay to two photons each, producing a large multiplicity of soft photons not unlike the 



13 



photon signal we have been discussing thus far. The remaining two thirds will consist of soft charged 
pions which will leave a large number of soft tracks in the tracking systems of the LHC experiments. 
We thus summarize that both hadronic and "photonic" contributions to the underlying event that 
come from new physics may be discovered by either or both of the following searches: a large 
multiplicity of soft tracks, and soft and unclustered energy deposited in the calorimeters. 

3.2 The Soft Track Signal 

The efficiency for discovering an anomalous underlying events by tracking depends highly on the 
tracking system, the algorithm used for tracking and the soft tracking efficiency in the LHC 
experiments. Here we will merely estimate the expected number of soft tracks from new physics 
and from backgrounds, and identify potential experimental challenges. 

What is the number of expected soft charged tracks in standard model underlying events? This 
has been studied extensively for various event generators in preparation for the initial tuning of 
MC generators during early running of the LHC. The numbers may vary among different event 
simulations, but it is useful to get ballpark estimate. Eventually these quantities will be measured 
and compared to signal (as describe in the strategy in section 1). For example, a study which 
is useful for our purpose [7] counts the number of soft charged particles in the central region in 
Drell-Yan di-muon events. For hard di-muon pairs, above a TeV or so, the number of charged 
particles per unit rj and 4> with |?7| < 1 is about dNch/drjdcf) ~ 0.7 — L7, depending on the Monte 
Carlo generator used. To get the total number of charged particles we multiply by An giving of 
order 10-20 tracks. At high luminosity one should add a significant number of tracks from pile-up 
of soft minimum bias events. At 4 x 10^"^ cm~^ s~^ the expected number of min bias events per 
bunch crossing is of order 10 (a cross section of ~ 100 mb at 0.1 mb~^ per bunch crossing), and 
more than double that at design luminosity. Each MB event is expected to produce 10-16 charged 
tracks in the central region [8]. In total, a conservative tally may produce as many as 200 particles 
per event in the |?7| < 1 region. 

How does our signal compare? The total number of photons is given by 



where f-y is the fraction of the energy emitted in photons. We estimated above that roughly 30% 
of the outgoing soft photons convert to electron-positron pairs in the tracking system. Taking a 
quirk mass of 500 GeV with A = 10 GeV and assuming to be 10% we get about 200 conversions 
in the tracking system (and thus about 400 tracks). Due to the geometric pattern it is reasonable 
to expect an order one fraction of these tracks to be in the central region. This is a reasonable 
starting point. The ability to claim that there is an excess of soft tracks from new physics will 
depend on how well the typical underlying event will be measured, the standard deviation from the 
average underlying event, and the exact size of the backgrounds, as well as the amount of energy 
emitted in photons versus hidden glueballs. In this search one may also search for the geometric 
pattern of the signal photons. This may be particularly useful since the charged particle density 
from backgrounds is expected to be flat in the r] — (p plane. 

We should stress however that the tracking signal that originates in photons may suffer from 




(11) 
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serious systematic issues. For example, it is not clear that a soft photon that converts somewhere 
in the middle tracker will be identified as two tracks. The conventional soft track algorithms may 
reasonably require some hits in the inner-most layers of the tracking system [7] , whereas a displaced 
photon conversion may "skip"' these inner layers. This would reduce the 30% estimate above to a 
few percent. This concern perhaps implies that counting tracks may not be the most efficient way 
to search for our signal. Instead one could hope to construct an observable that that gauges the 
amount of activity in the tracking system as a function of rj, (p, and perhaps radial depth r. The 
tally of charged tracks above should thus be taken as a rough indication for the feasibility of the 
search and its challenges. 

For colored quirks, on the other hand the prospects are more promising. Most of the energy is 
lost by soft pions, an order one fraction of which will be charged (say two thirds) and leave soft 
tracks. In this case these will be genuine tracks that originate from the interaction point, rather 
than converted photons. In this case the number of pions will be roughly 

CO 3 TT 3A2 ^ ^ 

For a quirk mass of 500 GeV and A ~ 10 GeV this gives of order a thousand charged tracks. This 
seams promising, given the expected backgrounds. Here too, however some systematic issues may 
be important, such as the tracking efficiency at such high multiplicities. Furthermore, investigations 
with a toy detector simulation (see next section) show that some of the charged pions may reach 
the EM calorimeters. These soft charged pions, which are highly curved, may "pollute" the pattern 
of soft energy in the Ecal which will be studied in the next section. 



4 Anomalous Underlying Events in the LHC Calorimeters 



In this section we will consider the feesability of observing anomalous underlying events using the 
distribution of soft energy in the electromagnetic calorimeters. We will simulate background and 
signal events. Due to the potential sensitivity to matter effects in the detector we simulate a "toy 
detector" using GEANT4 which reproduces the key features required for our analysis. Finally, we 
will analyze a sample of signal and background events and propose a method to distinguish them 
using a "CMB-like" multipole decomposition. We will propose two possible statistical analyses, 
one of which is gear towards a signal in soft photons and the other towards hadrons. 



4.1 The Signal and the Background 

For the purpose of our analysis we generated 500 signal events and 500 background events. 

Background: Our background sample should consist of the underlying part of hard VF+photon 
events. A significant portion of this background is expected to come from pile-up of minimum 
bias which are typical soft QCD events. Therefore, as our background sample we generated 500 
modified minimum bias events as follows. Min-bias events were generated using PYTHIA [19] 
tune A, including pile-up at a luminosity of 4 x 10^^ cm^^s"^ (0.1 mb~l per bunch crossing). 
Due to our limited ability to simulate the calorimeter, we modified the minimum bias by defining 
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all outgoing particles to be photons with a momentum that is identical to that of the original 
particle. A more in-depth analysis (i.e. one which employs the full detector simulation) would 
not make this modification. This modification, however, is very conservative in that it makes that 
background maximally similar to our calorimctric signal. In particular, in unmodified min-bias 
events soft charged particles will curve in the magnetic field and will not reach the EM calorimeter 
thus reducing the background in the Ecal. For this reason, it is important to keep in mind that a 
more careful analysis may show that the analysis of the upcoming section is underestimating the 
efficiencies for distinguishing the anomalous underlying event from the SM underlying event. 

Signal: As explained above, our signal is a large number of soft photons with an angular distribution 
similar to that depicted in Figure 1 . The details of the angular distribution and frequency spectrum 
in the center of mass frame are presented in Appendix A. We generated events according to this 
distribution. The total amount of energy in photons is {-s/S — 2mqi)f^ which is typically of order 1-2 
time f-yniqi. The fraction accounts for energy lost to invisible glueballs and is taken to be 1 or 0.1 
as discussed in Section 2.2 and s is generated according to the squirk production cross section. The 
photon distribution is further rotated and boosted longitudinally according to the differential cross 
section for quirk production. When generating these events we used the MRST parton distribution 
functions [27]. It is interesting to note that the longitudinal boost did not affect the pattern of 
photons significantly because the heavy squirkonium is produced near rest. Obviously, these soft 
photons will be generated on top of the standard underlying event at the LHC. Therefore, to each 
signal event we added it own background distribution of soft photons. To do this we added to the 
signal a modified minimum-bias event which was generated as described above for the background. 

4.2 Simulating a Toy Detector 

The soft nature of our signal makes it particularly sensitive to detector effects. Several simple 
detector simulations exist and are freely available - examples of these are ATLFAST [15] and PGS 
[16]. However, none of these correctly capture the low energy physics relevant to the processes we 
are considering. The crucial question is: given a particular distribution of photons, how much of 
the energy reaches the calorimeters and to what extent is the angular pattern recognizable? 

The back-of-the-envelope estimate of section 3.1 shows that indeed a significant amount of 
energy will reach the calorimeters. Understanding the effects of showering in matter and a magnetic 
field on the pattern, as well as a systematic comparison with backgrounds require simulating the 
events and the detector. To do this reliably one would need to know how much material each 
photon traverses on its way to the calorimeter and how it showers. This, of course, is what the full 
detector simulations do. However, for the purpose of our initial theoretical study we can ignore 
the details of the inner detector. Instead we chose to roughly reproduce the material budget of a 
"typical" detector at various distances. For this we constructed a hypothetical toy detector (we 
arbitrarily designed it ATLAS-like) in GEANT4 [18], a general package for modeling the passage 
of particles through matter, putting together a toy tracking system and calorimeter. 

A cross section of our toy detector is shown in Figure 4. We model only the barrel region of 
the detector, extending out to 0.6. We focus on the barrel region because the material budget 
before the calorimeter is significantly increased for the end caps. The inner tracker of our toy 
detector is represented by concentric cylinders of silicon - these represent the pixels and SCT of 
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Coil+Cryo 




Figure 4: A cross section of the toy detector used for our analysis. This simple detector (inspired by 
ATLAS) is sufficient to answer the simple questions - how much electromagnetic energy is expected 
to reach the various parts of the calorimeter? and is the angular pattern of soft energy visible? 
Five sample showers are displayed above, originating from five low energy photons. Evidently, four 
of these photons showered in the coil and one converted to e'^e~ in the tracking system. Only 
charged particles are shown. 

ATLAS, but for us the division into pixels and micro strips is ignored. Beyond the inner tracker 
is an outer tracker consisting of Argonne gas representing the ATLAS TRT. The total size and 
material budget of our tracking system is similar to that of ATLAS. Between the TRT and the 
calorimeter lie the coil and cryostat. This region is responsible for most of the energy loss between 
the interaction vertex and the calorimeter. Ignoring the detailed structure of these components, 
we modeled them by a single cylinder of aluminum with a radial width normalized to produce the 
correct material budget in front of the calorimeter. The calorimeter itself is modeled by a single 
layer of liquid argon cells divided along the (j) and z axes. Finally we placed the entire system 
in a 2T magnetic field. Our toy calorimeter is segmented into 20 bins in the z direction, ranging 
from -800cm to 800cm and 20 bins in the (j) direction. Our binning is quite arbitrary. The actual 
analysis of the diffuse energy deposition in the ATLAS calorimeter will involve the identification 
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of energy clusters using a Topocluster algorithm^ However, such an analysis is beyond the scope of 
this paper. 

As an example, in Figure 4 we also show the propagation of five low energy photons through the 
inner detector to the calorimeter. The multiple scattering of the particles, as well as the bending 
of charged particle trajectories are evident. Note that only one of the photons showered in the 
tracking system while the majority showered in the coil and cryogenic system. This is consistent 
with the relatively low material budget in the tracking system. 

While this toy detector does not, by any means, replicate the full architecture of the ATLAS 
detector, it does capture the most important points. First, the energy loss through pair production 
of primary photons and subsequent radiation and absorption of the secondary particles is correctly 
reproduced. We verified this with a uniform distribution of photons. In that simulation, the average 
energy loss per photon was in good agreement with Equation 10 as well as with an equivalent run 
of the full ATLAS detector simulation^. Second, the EM showers produced by the photons as 
they propagate through the detector are fully simulated by GEANT4, producing a realistic angular 
smearing of the energy deposition. 

4.3 Results 

In Figure 5 we show energy deposition in the calorimeter for three sample events. The first is a 

signal event for which all the bound state energy was emitted in photons. The second is a similar 
event but with only 10% of the energy in photons, and the third is a (modified) minimum bias event. 
The generation of these events is described in section 4.1. The antenna pattern, two back-to-back 
"doughnuts" is clearly visible in Figure 5(a), where all of the energy is emitted in photons. When 
only a tenth of the energy goes to photons the pattern is degraded to two back to back blobs and 
a quantitative statistical analysis, as described below, may be needed to label such an underlying 
event as an anomalous one. 

One could have worried that any angular pattern in the calorimeter would be smeared by 
magnetic field effects. However we already see that this will not be the case. This is because the 
soft photons may be divided into two rough groups - ones that convert in the tracking system 
and those that do not. The e+e^ pairs that arc produced by the early converting photons are 
highly curved by the magnetic field and thus barely deposit any energy in the calorimeters^. On 
the other hand, the photons that contribute most to the calorimetric signal are the ones that pass 
the tracking system without converting and are unaffected by the magnetic field. 

In order to claim an excess of anomalous underlying events within the signal region one would 
like to quantify the most atypical features of the signal events, perform cuts on the data and show 
a statistical excess compared to the expected background. One simple possibility is to cut on the 
total amount of unclustered energy in the barrel region, rj < 0.6. When all of the energy is emitted 

^This algorithm for collecting data from the calorimeters is different than the usual "sliding window" algorithm 
for identifying photons. In particular the topocluster has a much lower threshold, of order 100 MeV in the barrel 
region, which is set by noise [17]. 

'^We thank Elliott Cheu for this important cross check. 

®This is especially true at ATLAS where the coil is before the calorimeters. 
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Figure 5: Calorimeter energy deposition in the toy detector simulation. The distribution is shown 
for (a) bound state radiation with 100% of the energy released in photons, (b) bound state radiation 
with 10% of the energy in photons and (c) a minimum bias event. Brighter squares indicate a higher 
energy deposition in the cell, however, the scale itself is arbitrary for each figure separately. 



in photons, the average amount of energy deposited in our toy calorimeter is approximately 550 
GeV for a squirk mass of 500 GeV. For comparison, in the average (modified) min-bias event the 
average was below a 100 GeV. Our modification of the min bias events (see section 4.1), which was 
geared toward generating conservative backgrounds for pattern recognition (see below), may have 
increased the later number, but it may be taken as a ballpark figure. 

In this work we will focus on amore distinct "smoking gun" feature of our signal, the angular 
"anntena" pattern of soft energy. Identifying this pattern provides a unique data analysis challenge 
since most triggering and clustering algorithms are geared toward the identification of hard objects. 
A promising way to quantify the angular correlations of any function defined on a 2-sphere is to use 
a multipole decomposition, as was shown to be very effective in studies of the cosmic background 
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Figure 6: Multipolc expansion of calorimeter energy distributions. The multipole coefficients are 
shown for the bound state radiation (sohd) and the minimum bias event (dashed) from Figure 5. 
The normalization of the y axis is arbitrary. 

radiation^. Given a function f{9,(f)) The Z*'* multipole coefficient is 



The coefficient Ci receives its main contribution from fluctuations on angular scales of ir/l. Thus 
for the radiation pattern of Figure 5(a) we expect the I = 2 coefficient to dominate. This is 
exactly the effect that is visible in Figure 6, with higher multipole moments providing sub-dominant 
contributions to the expansion. The dashed line in Figure 6 represents the multipole decomposition 
for the minimum bias event shown in figure 5(c). The difference is quite evident. 

While the difference in the multipole decomposition between specific signal and minimum bias 

events is striking, it is not statistically significant on it's own. In fact, the minimum bias events 
exhibit a wide variation and many of them are atypical and could be mistaken for bound state 
radiation. In particular, it is known (and confirmed by our background event sample) that some 
underlying events will contain two broad back-to-back jets which may also create a peak at Z = 2, 
as well as higher even Vs. It is thus useful to quantify the expected peak at Z = 2 and look for an 
excess of signal above background. 

There are many observables which can be constructed from the first few multipoles. One 

^This analogy may be taken further, comparing the subtraction of the galactic plane to the removal of the end 
caps in our analysis. 




m 



where 
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possibility which we use for our statistical analysis is the following variable: 



sma = ■== , {16) 

where the averages are taken over the first five multipoles. If we consider these multipoles to be 
the components of a vector in a five-dimensional space, then sina is precisely the sine of the angle 
between this vector and the diagonal passing through the origin. An event in which all multipoles 
have the same magnitude will have sin a = whereas an event which is purely in a single multipole 
has sina = 1. Hence, this variable measures the degree to which the multipoles tend to differ from 
each other. Since the signal events are characterized by large differences between the I = 2 and the 
other multipoles, we expect sina to be larger compared to the minimum bias events which tend to 
be flatter. 

In Figure 7 we compare the sina values for three samples of 500 events. The first sample 
consists of signal events (overlayed on top of background UE's) , with all the energy of the bound 
state emitted in photons (solid), the second consists of similar events for which only 10% of the 
energy was emitted in photons (dashed) and the third are the background events (dotted) . Imposing 
a cut of sina > 0.7 we find efficiencies for the three event samples as listed in Table 1. The variable 
we used above does not single out any particular multi-pole. Other possible observables that focus 
on Z = 2 might be used, though we did not find that these did much better in our particular case. 



200 




0.4 0.6 

sina 



Figure 7: sina for samples of 500 events: signal events with %100 of the energy emitted in photons 
(solid), 10% of the energy in photons (dashed) and minimum bias events (dot-dashed). 
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% of energy 


Efficiency for 


in photons 


sin a > 0.7 


100 


0.74 


10 


0.17 


Background (0) 


0.02 



Table 1: Efficienciey for passing a cut of sin a > 0.7 for the three event samples of 100% of the 
energy in photons, 10% of the energy in photons, and standard model backgroud. The variable sin a 
measures to what extewnt the angulard distribution of soft energy in the calorimeter is dominated 
by a single multipole. 

4.4 Analyses I: Discovering Anomalous Underlying Events 

Given the efficiencies above, one may analyze the data by several approaches depending on the 

efficiency and confidence level of the original hard search and on the amount of energy emitted in 
soft radiation with each new physics event. Wc will begin with a scenario that is quite likely in our 
case of uncolored squirks annihilating to W'^. This hard final state is quite distinct and the signal 
to background ratio may well be of order one, as demonstrated in [13]. However, even when new 
physics is discovered we would like to extract as much information about it as possible. In other 
words, would like to test whether the discovered signal is a "normal" resonance in W'y or perhaps 
a quirky bound state. We would like to test whether the new signal is correlated with anomalous 
underlying events. 

To test this hypothesis we may exploit the fact that the underlying event is believed to be 
quite universal, particularly for similar final states. Following the strategy outlined in section 1 
we can thus collect a "control" data set of standard underlying events by considering similar final 
states in different kinematic regions, or completely different final states (e.g. di-muon events). In 
fact, the characteristic properties of underlying events, such as charged particle density and px 
distributions have been shown to be independent of a/I in jet events as well as Drell-Yan (beyond 
a certain energy) [7] . This implies that one can collect a sizable pool of standard underlying events 
which arc not "contaminated" by contributions from new physics such as quirks. This large control 
sample can then be compared with the underlying events in the signal region. 

Due to the large number of events, such an analysis may be effective in discovering quirky 
dynamics even under the more pessimistic assumption of 10% of the excitation energy being released 
in photons. Given enough statistics the the difference between the 10% and BG distributions of 
figure 7 may be established. 

To see this lets do a gedanken-analysis with 100 fb^^ of data and a simple two-bin chi-squared 
test. In line with the analysis of [13] we will assume the quirk production cross section is of order 10 
fb and the branching ratio to leptons times efficiency for the search is of order 0.1, giving roughly 
100 signal events. These numbers arc consistent with a squirk mass of 500 GeV, which was used in 
our simulations. With these efficiencies the signal to background ratio is of order 1, providing clear 
evidence for new physics well beyond 5 sigma. Thus the "signal sample" consists of 200 events, a 
half of which are true signal. The multipole variable sin a in signal and background events follows 
the distributions shown in figure 7 (of course, a more detailed analysis of the UE background and 
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the detector may yield different results). As a simple example consider just two bins of sin a bellow 
and above 0.7. Our simulation showed that the fractions of events with sin a > 0.7 are 0.02 and 
0.17 for background and signal events (10% of energy in photons) respectively. On average, about 
20 of the 200 events will fall in the since > 0.7 bin. 

This may be compared with a large control sample consisting of pure background underlying 
events. In a control sample of 3000 events only 60 events will fall in the sin a > 0.7 bin on average. 
If the measured frequencies landed close to these averages a chi-squared test of independence gives 
a very "bad fit" of ^ 40 for a single degree of freedom. Such a measurement would support that 
the two samples are drawn from two different distributions with a high confidence level. 

The ability to consider larger control samples is important to stabilize such a result against 
statistical fluctuations. For example, consider the case in which we are unlucky and the actual 
measured numbers of events in the sin a > 0.7 bin fluctuate by 2 sigma down and up for both 
the signal and the control samples respectively. Even in this unfavorable case the measured signal 
and control samples are statistically independent with ^ ^ ^ single degree of freedom. This 
demonstrates that in order to carry out this program and large control sample of standard model 
underlying events should be collected and studied. 

4.5 Analysis II: Improving the Confidence of Discovery of New Physics 

We now consider a different situation (which is not necessarily applicable in our case of uncolored 
squirks going to W'y). Lets assume that an order one fraction of the excitation energy (of order 0.5 
TeV) is released in photons. However, we can also suppose that the confidence level of discovery 
in the standard hard analysis is poor, bellow 3 sigma. In this case we can exploit the distinct 
soft signal (as seen in the 100% curve of figure 7) to improve the confidence level of discovering 
new physics. Instead of searching for the hard final state X, we should perhaps search for X+an 
anomalous underlying event. Such an analysis may fit well with the case of colored quirks in which 
the amount of soft energy is perhaps larger (since the competition with invisible glueballs is less of 
an issue) and the hard final state may be in jets, which suffers from a larger background than the 
distinct W^. As we pointed out earlier the systcmatics of soft hadrons may be somewhat different 
than those of and should be considered separately. 

To demonstrate this lets again assume the efficiencies derived from figure 7. We can now use 
the atypical shape of the underlying event as a discriminant between signal and background. In 
Imposing a cut of sin a > 0.7 the background will be efficiently reduced to 2% of that in the 
original search whereas the signal is reduced by ~ 0.74. Therefore, in the more optimistic scenario 
of 100% of the energy emitted in photons, the original S/-\/B may be improved by a factor of 
0.74/V0.02 ~ 5.2. In this case, our underlying event analysis may potentially take a 1-2 sigma 
bump into a discovery beyond 5 sigma! The less optimistic case of 10% of the energy in photons 
the improvement factor is a more modest 1.2. 
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5 Discussion 



In this work we have focused on a specific example of new physics, uncolored squirks, which lead 
to a diffuse soft signal at the LHC. As we mentioned throughout the paper colored quirks may 

also be promising candidates for generating soft signals in the spirit discussed here. Are there 
any other cases? The Hidden Valley scenario [1] has been proposed recently as an interesting 
possibility for novel LHC signals. The general idea involves the production of new particles which 
are charged under a new strong force. These particles will shower and hadronize according to the 
strong dynamics of the new sector. Such a showering process will typically distribute the high 
energy of the produced particles into many soft valley hadrons, some of which will decay to soft 
SM jets, leptons or photons. One promising way to produce a diffuse soft signal in the spirit of the 
ones we are studying is to consider an unparticle scenario [23] in which the conformal symmetry 
in the hidden sector is broken at some low scale^'^. The strong conformal dynamics will produce 
very broad "jets" of valley particles (the part of the spectrum just above the gap, which is likely 
to contain resonances) and thus broad distributions of soft and diffuse SM decay products. Indeed, 
"multi-unparticle " production may be enhanced in these scenarios [25] . 

One particularly interesting possibility is that the hidden sector couples to the SM through the 
"Higgs portal" [26, 2]. For example, the Lagrangian may contain an operator 

C D \H\^Ov (14) 

where Oy is an operator in the new strong sector. If the hidden sector indeed contains strong 
and nearly conformal dynamics, the most likely decay channel of the Higgs may be to broad and 
soft "jets" of valley hadrons. The decay of the light valley particles to standard model states will 
produce a soft modification to underlying events at the LHC. A possible search would be for a single 
boosted Z plus an anomalous underlying event from the Higgs decay. The multiploe expansion of 
soft energy in this case would presumably exhibit a peak at Z = 1 with subleading peaks at higher 
odd Vs. 

Another interesting possibility which has been proposed recently is that the Higgs itself is part 
of the unparticle sector [28] . In this scenario the Higgs is part of a strongly coupled sector which 
also contains very hght degrees of freedom. These light degrees of freedom may couple back to the 
standard model, e.g. by integrating out the heavy Higgs. It is thus plausible that as a Higgs is 
produced a shower of light degrees of freedom is emitted which consequently decay back to soft SM 
particles. Such an event will contain the hard decay products of a Higgs in addition to a modified 
underlying event. In these scenarios, the strategies outlined in this work may ultimately shed new 
light on the nature of the Higgs. 

In conclusion, in this work we have proposed a new signal of new physics at the LHC. Searches 
for new and heavy physics has rightfully focused on hard physics object. We have demonstrated 
that new physics may also be discovered in the soft, or underlying component of LHC events. 
Discovery of an anomalous underlying event may be a unique way to discover the non-perturbative 
nature of new physics that is otherwise unattainable. We demonstrated our strategy in a simple 
model of quirks (or squirks) which is inspired by folded supersymmetry, a model for addressing the 

^"This scenario may thus be called a hidden valley [24]. 
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electroweak hierarchy problem. We have discussed some of the strategies to search for anomalous 
underlying events, both in the tracking and the calorimetry of the LHC experiments. 

A general lesson to be drawn from this work is that underlying event studies should accompany 
some of the standard searches of new physics. Such an underlying event study may either enhance 

the confidence of discovering new physics, or teach us new information about new physics. The 
hard searches which may benefit from an accompanying UE analysis include searches for heavy 
resonances, (di-boson, di-lepton or di-jet) as well as searches involving missing transverse energy. 
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A Classical trajectory and Photon Spectrum 

To generate the photon spectrum in our simulation we modeled the radiative decay of squirkonium 

by the classical radiation of accelerating charges. The charged squirkonium is modeled as two 
particles of equal masses and of charges ±2/3 and ±1/3 with a constant force acting between them. 
The classical equation of motion is 

= A^ (15) 



dt 

in the center of mass frame. Since the squirks are produced semi relativistically we use relativistic 
momentum p = mj/d. Assuming the motion is along the x direction the solution is a periodic 
trajectory with a period of 

^=^^2 (16) 

with the trajectory for <t < T/2 described by 

E - Ve"^ - 2AV-E2 - Am^t + 4AH^ , , 

= ^ ^ 

where E is the total energy of the system, initially set to V^. The classical photon spectrum 
may be evaluated by Fourier decomposing the retarded field far away from the oscillating charges. 
Because of the periodic nature of the motion the field decomposes into a discrete Fourier series 
with frequencies = nwo and ujq = 27r/T. Following [22] , the power radiated per solid angle in 
the nth frequency mode is 



27r/a)o ^ 

" inijJo{t+Ri(t)) 



(18) 



where the i sum is over the two charges, qi is the charge of the particle i and Ri{t) is the 
instantaneous distance of a far away point from that particle. Assuming the linear motion of 
equation (17), the Fourier integral over the path of the particle may easily be done numerically. 
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To simulate the radiative decay of a quirkonium system produced with a center of mass energy 
of \f§ we simply generate a photon according to the spectrum of equation (18) using a trajectory 
of that energy. We then subtract the photon's energy from the oscillating system and generate the 
next photon with the reduced energy, and so on. As discussed in Section 2, one should assume 
some of the energy is radiated in the form of (hidden) glueballs. For example, as a benchmark 
we also consider the case where only 10% of the energy is radiated in photons. To generate these 
events we follow the same procedure, but only considering every tenth photon as visible, throwing 
90% of the photons away. 
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